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Section 1 
INTRODUCTION 


The objectives of this study were to provide both qualitative and 
quantitative comparisons of several radiometric processing procedures 
applied to MSS data, to assess the adequacy of present MSS calibration 
procedures, and to assess the impact of a modified MSS radiometric pro- 
cessing procedure on the throughput of the Master Data Processor. 

The radiometric calibration currently performed by NASA on image data 
received from the Landsat Multispectral Scanner (MSS) is based on a 
linear radiometric transformation derived dynamically from six samples 
of the calibration wedge data provided by the MSS for each image line. 
This radiometric calibration procedure leaves a residual periodic 
variation in image radiometry which is termed striping. This radio- 
metric variation, in addition to its impact on the visual quality of 
processed imagery, introduces a systematic radiometric error which may 
compromise the validity of results obtained by machine processing for 
information extraction. 

Over the years, various techniques have been proposed to compensate MSS 
digital data for striping. These proposals have concentrated on tech- 
niques which characterize the response of individual MSS detectors by 
examining the image data reported by each, and equalize this response 
across the detectors of a spectral band. These techniques have been 
shown to provide considerable reduction in visible striping. They fail, 
however, to address the questions of why the calibration procedures 
employed are not effective in eliminating striping originally, and 
whether the calibration procedures can be modified to eliminate striping. 
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These questions could only be answered by studying the calibration 
procedures and the data available to perform the calibration- It was 
anticipated that investigation of the temporal stability of the complete 
set of calibration wedge data , together with interdetector comparisons 
of this data, would provide the information necessary to formulate a 
more effective MSS calibration procedure, thus eliminating these inter- 
detector variations as part of the calibration process. 
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Section 2 

PROCESSING DESCRIPTION 


The study was designed to empirically investigate both destriping pro- 
cedures and the MSS calibration procedure, evaluating their effectiveness 
by both visual and quantitative procedures . It also considered the 
characteristics of MSS calibration wedge data in order to elucidate the 
temporal and interdetector effects which exist in this data. 

The study consisted of the following tasks : 

a. Preprocessing of scene CCT to produce band separated images 
with associated calibration data 

b. Destriping processing and MSS calibration processing 

c. Visual evaluation of processing results 

d. Power spectrum analysis of processed image data 

e. Detector-specific image histogram generation 

f. Multispectral classification of agricultural areas in processed 
scene data 

g. Evaluation of MSS Calibration Data 

h. Evaluation of the impact on MDP throughput of the selected 
destriping procedure. 
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One Landsat scene> provided in both calibrated and uncalibrated form, 
was used in this study. This MSS scene, E-2183-16433, contains water, 
rural crop areas, and areas of high and low radiance. Both the radio- 
inetrically calibrated and radiometrically uncalibrated versions of this 
scene were processed and evaluated as described below. 


2 . 1 PREPROCESSING 

Each of the two versions (calibrated and uncalibrated) of the scene were 
reformatted to produce complete images of each of the four MSS bands, 
with the calibration data from the CCT provided for each MSS mirror 
sweep. These images were radiometrically adjusted to match the charac- 
teristics of the IBM Drum Scanner/Plotter, and recorded on photographic 
film, to permit inspection for data quality. 


2.2 RADIOMETRIC PROCESSING 

The radiometric processing techniques investigated in this study are 
outlined below. 


2.2.1 IBM's Destriping Procedures 

Type al. Mean and standard deviation equalization based on the 

data from a single MSS mirror sweep (i.e., six contiguous 
lines of image date) , with the resulting detector unique 
gain and bias radiometric compensations applied to the 
data from the subsequent mirror sweep. 
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Type a2. Mean and standard deviation histogram equalization based 
on the data from the first 60 mirror sweeps of the image 
(approximately 15% of the image data in an MSS frame) , 
with the resulting detector unique gain and bias radio- 
metric compensations applied to the data for the entire 
image . 

Tjrpe a3. Mean and standard deviation histogram equalization based 
on the data from the first 195 mirror sweeps of the image 
(50% of the image data) , with the resulting detector 
unique gain and bias radiometric compensations applied to 
the data for the entire image. 

Details of this processing are presented in Appendix A. 

2.2.2 Modified MSS Calibration Procedures 

Type bl. The smoothed gain and offset values developed for the 

lines of the initial 98 mirror sweeps of a spectral image 
(approximately 25% of the image data) were averaged for 
each of the six detectors in a spectral band, and these 
six averaged gain and bias sets were used as the radio- 
metric correction for the entire image. This smoothed 
gain and offset calibration employed the six calibration 
wedge values present in the CCT data for the scene. 

Type b2. The smoothed gain and offset values developed for the 

lines of the entire spectral image (100%) were averaged 
for each of the six MSS detectors of that band, and these 
six averaged gain and bias sets were used as the radiometric 
correction for the entire image. 

The algorithm used to calculate the smoothed gain and offset is presented 

in Appendix B. 
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Section 3 

EVALUATIONS OF PROCESSED IMAGES 


The comparative evaluations which were performed in this study of 
radiometric processing procedures are shown, in Figure 3-1, together with 
the data processing required to generate products for comparison. These 
results of evaluations are discussed in the following subsections. All 
the evaluations except raultispectral classification were performed on 
the scene data from MSS band 1 (0,5 to 0.6 micrometers), which exhibited 
the greatest visually apparent striping in photographic recordings of 
the unprocessed, uncalibrated scene data. The multispectral classification 
evaluation employed data from all four MSS bands. 
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3.1 EVALUATION BY DETECTOR-SPECIFIC HISTOGRAMS 


A comparative evaluation among images of MSS band 1 was made using 
radiometrically uncalibrated data processed by the type al, a2, and a3 
destriping procedures and the type hi and b2 MSS calibration procedures , 
and radiometrically calibrated data processed by the type al destriping 
procedure. The evaluation was made on the basis of histograms of the 
data reported from each of the six MSS detectors in each of the six 
radiometrically processed spectral images. Detector-specific and aggregate 
histograms of this data were generated, and the consistency of each 
individual detector histogram with the corresponding aggregate histogram 
was considered. These histograms are presented in Figures 3.1-1 through 
3.1-6. 


A conventional statistical approach to evaluating the goodness of fit 
between the aggregate histograms and their associated detector-specific 
histograms was used to test the hypothesis that the frequency distribution 
of pixel values for an individual detector is one sixth of the frequency 
distribution of pixel values over all detectors, using the statistic: 


Where 

N is the number of occurrences of pixel value i for all 
1 

detectors . 


k 

= E 

i«l 


(N^/6) 


N is the number of occurrences of pixel value i for 
Ji 

detector j , and 


All of the calculated X. values for the histogram sets, as given in 

J 2 

Table 3.1-1, greatly exceed the critical X value at the 0.005 level of 
significance (i.e., there is less than a 0.5% probability that a random 
set of pixel values, taken from a population whose frequency distribution 
is given by the aggregate histogram, would show such poor agreement with 
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the aggregate histogram) , forcing the conclusion that none of the pro- 
cessing techniques has achieved a completely satisfactory equalisation 
across the detectors. 

2 

The calculated values do, however, provide a summary measure of the 
deviations of the distributions in each of the detector-specific histo- 
grams from the distribution in the corresponding aggregate histograms. 

2 

By summing the X. values for the six detectors of a spectral image, an 
overall measure of these deviations is obtained which can be used to 
rank the performance of the various radiometric processing techniques. 
On this basis the best equalization of the detector histograms was 
obtained with the type a2 processing of the uncalibrated scene data. 
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Table 3.1-1. Chi-squared Values for Comparison of Detector-Specific Histograms 


Radiometric 

Processing 

Degrees 

of 

Freedom 

in 

Values 

J 

Sum 

0.005 

3 = 1 

II 

3 = 3 

3 = 4 

3 = 5 

j = 6 

Uncalibrated Scene Data 
after Type al Processing 

10 

254587 

1316226 

375530 

9126 

39104 

8929 

2,003,503 

25.2 

Uncalibrated Scene Data 
after Type a2 Processing 

10 

36392 

50097 

65716 

47192 

127396 

52941 

379,733 

25.2 

Uncalibrated Scene Data 
after Type a3 Processing 

10 

35113 

178514 

41636 

80355 

102538 

80897 

519,054 

25.2 

Uncalibrated Scene Data 
after Type bl Processing 

17 

487846 

69935 

174762 

67419 

92641 

63042 

955,646 

35.7 

Uncalibrated Scene Data 
after Type b2 Processing 

17 

488333 

69918 

175037 

67403 

92638 

63020 

956,348 

35.7 

Calibrated Scene Data 
after Type al Processing 

19 

130507 

104505 

259655 

341799 

417700 

421118 

1,675,284 

38.6 
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Figure 3.1-3. Detector-specific Histograms for Uncalibrated Scene Data 
from MSS Band 1 after Type a3 Radiometric Processing 
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Figure 3.1-4. 


Detector-specific Histograms for Dncalibrated Scene Data 
from MSS Band 1 after Type bl Radiometric Processing 
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Figure 3.1-5. 


Detector-specific Histograms for Uncalibrated Scene Data 
from MSS Band 1 after Type b2 Radiometric Processing 
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Figure 3.1-6. 


Detector-specific Histograms for Calibrated Scene Data 
from MSS Band 1 after Type al Radiometric Processing 
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3.2 EVALUATION BY VISUAL INSPECTION 


A comparative evaluation among Images of MSS Band 1 was made using 
radioraetrically uncalibrated data processed by the type al, a2, and a3 
destriping procedures and the type bl and b2 MSS calibration procedures. 
The method of evaluation was visual inspection of film recordings of the 
radiometrically processed images . The images were ranked on the basis 
of the visually apparent striping which they exhibited. 

Photographic prints of the scene data used in this visual evaluation are 
presented as Plates I through 6. The actual evaluation was made from 
the negatives used to produce these prints, in order to avoid any varia- 
tions which might be introduced in producing prints from the negatives , 
and to provide the best possible presentation of the data during evalua- 
tion. Negatives of image data for the five radiometrically processed 
versions of the scene for MSS band 1 (0.5 to 0.6 micrometers), together 
with a negative of the unprocessed, uncalibrated data, were arranged 
randomly in a two-wide by three-high array and uniformly illuminated 
from behind. Five experienced members of IBM’s digital image processing 
group were then asked to rank the six images on the basis of the visually 
apparent striping. 

The rankings which were obtained are presented in Table 3.2-1, together 
with an overall ranking obtained by summing the individual rankings. On 
the basis of this overall ranking, the tj^pe a2 processing of uncalibrated 
scene data is clearly preferred. 


3-12 



Table 3.2-1. Results of Visual Evaluation of Radiometrically Processed Scene Data 
for Landsat MSS Scene 2183 - 16433, MSS Band 1 


Image 

Number 

Image 

Description 

Rank; 

ngs Obtained* 

Overall 

Ranking 

Case 

1 

Case 

2 

Case 

3 

Case 

4 

Case 

5 

1 

Uncalibrated scene data, 
type bl radiometric processing 

1 

2 

4 

3 

3 

13 

2 

Uncalibrated scene data, 
type al radiometric processing 

1 

5 

5 

2 

1 

14 

3 

Uncalibrated scene data, 
type b2 radiometric processing 

1 

1 

1 

5 

5 

13 

4 

Uncalibrated scene data, 
type a2 radiometric processing 

1 

3 

2 

1 

2 

9 

5 

Uncalibrated scene data, 
type a3 radiometric processing 

1 

4 

3 

4 

4 

16 

6 

Uncalibrated, unprocessed 
scene data 

6 

6 

6 

6 

6 

30 


* 1 = best 
6 = worst 


















3.3 EVALUATION BY POWER SPECTRA ANALYSIS 


A comparative evaluation among images of one spectral band was made 
using radiometrically uncalibrated data processed by the type al, a2, 
and a3 destriping procedures. The method of evaluation was power spectra 
comparison. 

Power spectra for two 128-line by 50 sample subimages were generated for 
each of these radiometrically processed images. The subimages were 
chosen to provide power spectra of a low frequency and a high frequency' 
regron. Power spectra were taken only in the "along track" direction of 
the image, since the linear gain and offset radiometric transformations 
employed by the destriping and calibration procedures have known effects 
on the power spectra in the along scan direction. The specific subimages 
are identified in Plate 7, which is an overlay for Plates 1 through 6. 

The power spectra obtained for this evaluation for regions in the 
uncalibrated scene data and in the three radiometrically processed 
scenes are presented in Figures 3.3-1 through 3.3-12. In these figures, 
the abscissa is the unnormalized squared value of the discrete Fourier 
coefficient, while the ordinate is the index of the coefficient, and 
ranges from 0 to the index corresponding to the folding frequency, which 
is 64 for the 128-line transforms used in this evaluation. 

In engineering terras, the spacing between the Fourier components is: 

■ (128 lines) (80 meters) ' 9.766 X lO’^ cycles/ncter 

while the frequency range of the plotted spectra is: 

= 64(AF) = 6.25 X 10 ^ cycles/meter. 
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Components of the power spectra which correspond to striping would be 
expected to occur at intervals of: 


■ 6 (SO'LtlS) ■ cycles/meter. 

corresponding to index values i = n [(1/480)]^/^ [(128) (80)] where n is 
an integer. For the frequency range of the plotted spectra, i = 21.3, 

42.7, and 64. 

In the power spectra of the high frequency region presented an Figures 3.3-1 
through 3.3-4, there is only minimal evidence in the uncalibrated scene 
data of greater power concentration at those frequencies which would 
correspond to striping than at other points in the spectra, and there is 
less evidence of this effect for the processed data. The peak in spectral 
power which is seen at the folding frequency, in the absence of lower 
frequency striping-related harmonics, cannot justifiably be attributed 
to striping, since power at this frequency alone corresponds to a two~ 
line periodicity rather than a six-line periodicity. In the power 
spectra of the low frequency region presented in Figures 3.3-5 through 
3.3-8, the expected concentration of power at the frequencies corresponding 
to striping is evident, in the unprocessed data, as is the suppression 
of power at these frequencies which results from the three different 
types of radiometric processing. 

A less complex power spectrum comparison of the striping in these images 
is provided by power spectra which were generated for a smaller (30 
sample by 64 line) low frequency region contained entirely within the 
originally defined low frequency region and containing only image samples 
of the Missouri River. These spectra are presented as Figures 3.3-9 
through 3.3-12. The simple power spectrum of this uniform- level subimage 
provides a background against which the power concentrations at the 
striping frequencies are much more apparent than they are for the power 
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spectra of Figures 3.3-1 through 3.3-8. For this smaller region, the 
Fourier coefficient indices calculated for the striping frequencies are 
10-7, 21.3, and 32.0, and power concentrations are quite apparent near 
these values in Figure 3.3-9, the power spectrum obtained for the uncali- 
brated scent data. 

A quantitative characterization of the power present at the striping 
frequencies for this small, low frequency subimage in the uncalibrated 
data and in the three processed data sets is presented in Table 3.3-1. 

In this table, the power levels at the striping frequencies are expressed 
in decibels above the average power level for the particular spectrum in 
order to normalize the data for this average, which differs among the 
source data sets. From this table, it is apparent that the type a3 
radiometric processing is most effective in suppressing the power at the 
fundamental striping frequency, and the type a2 radiometric processing 
IS marginally better than type a3 at suppressing the first harmonic. 

Both of these techniques considerably enhance the relative power present 
at the folding frequency. An exact explanation for this effect is 
unknown, but the detector-specific histograms presented in Figures 3.1-1, 
3.1-2, and 3.1-3 suggest that the "forbidden" pixel values produced by 
these two radiometric processing techniques, which are different for 
each detector, result in an effective increase in the quantization noise 
present in the image. Since such an increase would be detector-related, 
it could appear in the power spectriim at the highest frequency, the 
folding frequency. 

Overall, the type a3 radiometric processing provides the most effective 
combined suppression of the fundamental and first harmonics of striping, 
and is the most effective technique on this basis . 
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Figure 3.3-1. Power Spectrum for High Frequency Region in Uncalibrated Scene Data 
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Figure 3.3-2. Power Spectrum for High Frequency Region in Uncalibrated 
Scene Data after Type al Processing 
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Figure 3.3-3. Power Spectrum for High Frequency Region in Uncalibrated 
Scene Data after Type a2 Processing 
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Figure 3.3-5, Power Spectrum for Low Frequency Region 
in Uncalibrated Scene Data 
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Figure 3.3-6. 


Power Spectrum for Low Frequency Region in Uncalibrated 
Scene Data after Type al Processing 
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Figure 3.3-7. Power Spectrum for Low Frequenc 7 Region in, Uncalibrated 
Scene Data after Type a2 Processing 
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Figure 3.3-8. Power Specturm for Low Frequency Region in Uncalibrated 
Scene Data after Type a3 Processing 
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Figure 3.3-9. Power Spectrum for Smaller Low Frequency 
Region in Uncalibrated Scene Data 
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Figure 3.3-10. 


Power Spectrum for Smaller Low Frequency Region in Uncalibrated 
Scene Data after Type al Processing 
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Figure 3.3-11. 


Power Spectrum for Smaller Low Frequency Region in 
Scene Data after Type a2 Processing 
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Figure 3.3-12. 


Power Spectrum for Smaller Low Frequency Region in Uncalibrated 
Scene Data after Type a3 Processing 
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Table 3.3-1, Relative (to Average Power) Power Spectrum Values for Restricted 

Low Frequency Image Region at Frequencies Corresponding to Striping 



Fundamental 

f = 2. 15 X 10 cycles/meter 

First Harmonic 

I =321 

f = 4.10 X 10 cycles/meter 

Second Harmonic 
i °-32 

f = 6.25 X 10 cycles/meter 

Uncallbrated Scene Data 

4.04 db 

5.78 db 

4.83 db 

Uncalibtated Scene Data after 
Type al Radiometric Processing 

2.42 db 

1.80 db 

5.87 db 

Uncalibrated Scene Data after 
Type a2 Radiometric Processing 

0.08 db 

-0.08 db 

10.17 db 

Uncallbrated Scene Data after 
Type a3 Radiometric Processing 

0.02 db 

-0.09 db 

10.25 db 




3.4 EVALUATION BY MULTISPECTRAL CLASSIFICATION 


A comparative evaluation by multispectral classification was performed 
between two multispectral image data sets, one produced by application 
of the type bl MSS calibration procedures to the four bands of uncali- 
brated MSS data and the other produced by application of type al destriping 
procedure to the four bands of calibrated MSS data. 

Each four-band data set was evaluated using a parallelpiped (limit) 
classifier. For each data set, a training field for each of two classes 
was selected based on the ground truth data. These training fields were 
used to determine the class limits. A test field for each of these 
classes was identified, and the four-vectors within each of these fields 
were classified against the class limits for the class. 

This procedure was then repeated, exchanging the definition of training 
and testing fxelds. 

Plate 7 IS an overlay for Plate 5 which identifies the location in the 
image data sets of the Hand County LACIE Intensive Study Site 2, which 
contained the fields employed in this evaluation. Figure 3.4-1 is a 
diagram of the field structure in this site, and the specific fields 
employed (5, 50, 204, and 270) are identified on this diagram. 

The results obtained in these classification experiments are presented 
in Table 3.4-1. In all cases, the Type al radiometric processing on 
calibrated scene data has produced image data for which a higher percentage 
of the field pixels are correctly classified than has the Type bl radio- 
metric processing on uncalibrated data. However, if the evaluation is 
based on the volume of the four space parallelpiped which characterizes 
the classes, a preferred radiometric processing cannot be identified, 
since for spring wheat this volume is smaller for the Type al processing 
of calibrated data, while for oats this volume is smaller for the Type 
bl processing of uncalibrated data. 
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Because of conflicting conclusions which resulted from this experiment, 
and because of the lack of spearation of the class regions in four-space 
which is evident in Table 3.4-1, no reliable conclusions can be drawn 
for this evaluation - 
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Table 3.4-1. .Results of Evaluation by Multispectral Classification 


Class 

Field 

1 

ID 

MSS 

Band 

Class 

4-Space 

Range 

(Counts) 

Class 

Mean 

Vector 

(Counts) 

No of 
Pixels 

Field 

2 

ID 

MSS 

Band 

Class 

4-Space 

Range 

(Counts) 

Class 

Mean 

Vector 

(Counts) 

No of 
Pixels 

Radiometric 
Processing, 
Source Data 

Percent of Pixels Correctly Classified 
Number of Pixels Correctly Classified 

Freld 1 * Trarnmg 
Freld 2 = Test 

Field 2 = Training 
Field 1 = Test 

Mature 

Spring 

Wheat 

204 

1 

e 

23 5 

42 

270 

1 

11 

26 0 

30 

Type a 1 

Calibrated 

Data 

13'= 

4 Pixels 

71'-. 

30 Pixels 

2 

11 

29 0 

2 

34 

38 5 


3 

12 

53 5 

3 

18 

60 5 

4 

7 

28 0 

4 

9 

30 0 

4-S( 

aace Volume 5544 

4-Space Volum 

60588 

Mature 

Spring 

Wheat 

204 

1 

7 

24 0 

42 

270 

1 

14 

28 5 

30 

Type b 1 

Uncalibrated 

Data 

3 3% 

1 Pixel 

50% 

21 Pixels ! 

j 

i 

2 

11 

29 0 

2 

32 

44 5 

3 

11 

55 0 

3 

19 

62 0 

4 

7 

28 0 

4 

11 

31.0 

4-S( 

)3ce Volume 5929 

4-Space Volum 

e 93622 

Mature 

Oats 

5 

1 

13 

27 0 

37 

50 

1 

11 

26 0 

24 

Type a 1 

Calibrated 

Data 

83 3% 

20 Pixels 

67 6% 

25 Pixels 

2 

27 

37 0 

2 

20 

37 5 

3 

19 

62 0 

3 

23 

560 

4 

6 

30 5 

4 

11 

28 0 

4-Sp 

ace Volume 

40014 

4-Space Volume 55660 

Mature 

Oats 

5 

1 

8 

29 5 

37 

50 

1 

7 

29 0 

24 

Type b 1 

Uncalibrated 

Data 

54 2% 

13 Pixels 

541% 

20 Pixels 

2 

14 

44 5 

2 

10 

42 5 


3 

12 

65 5 

3 ' 

18 

60 5 


4 

6 

31 5 

4 

10 

28 5 


4-Spate Volume 8064 

4-Space Volume 

12600 
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Figure 3.4-1. 


Field Structure of the Hand County LACIE 
Intensive Study Site 
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3.5 SUMMARY OF COMPARATIVE EVALUATION RESULTS 


The results obtained in the comparative evaluations of the radioraetrically 
processed image data are summarized in Table 3.5-1, From this table, it 
js apparent that the Type a2 radiometric processing of uncalibrated 
scene data has provided the most satisfactory compensation for striping 
on an overall basis. Scene data produced by this technique was ranked 
first when evaluated both visually and by means of detector-specific 
histograms, and was ranked second when evaluated by power spectrum 
analysis. 

Since the multi. spectral classification evaluation was performed for only 
the type bl processing of uncalibrated data and the type al processing 
of calibrated data, the inconclusive results obtained from this evaluation 
do not compromise the identification of the type a2 processing of uncali- 
brated scene data as the best overall procedure. Had this evaluation 
identified a preferred processing technique, one could still not conclude 
that this processing technique was superior to the type a2 processing of 
uncalibrated scene data. 
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Table 3.5-1. Summary of Evaluation Results 
(1 = best) 


Lvaluation Method 

Uncalibrated Scene Data 
Processing Type 

Calibrated Scene Data 
Processing Type 
al 

al 

a2 

a3 

bl 

b2 

Detector Specific Histograms 

6 

a 

2 

3 

B 

5 

Power Spectrum 

3 

2 

1 

B 

jjBjlllllll 

- 

Multlspectrai Classification 

Inc 

oncli 

sive 

resul 

ts 


Visual Evaluation 

4 

i 

5 

2 

2 
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Section 4 


EVALUATION OF MSS CALIBRATION DATA 


The calibration wedge data for one MSS band was investigated to charac- 
terize its temporal stability, by detector, and the relationship between 
the calibration data for the various detectors. For this investigation, 
the full set of MSS band 1 calibration wedge data corresponding to the 
first 210 mirror sweeps of MSS scene E-2183-16433 was employed. 


4-1 



4.1 TEMPORAL VARIATION OF REPORTED CALIBRATION WEDGE VALUES 


For each detector of MSS band 1, plots of the mean and root-mean-squared 
variation of 50 calibration wedge samples were generated. These plots 
statistically summarized the temporal variation of a larger number of 
these values than the six points employed for calibration. These plots 
are presented as Figures 4.1-1 through 4.1-6. On these plots, the 
abcissa is the decompressed calibration sample value, and the ordinate 
IS the sample number, counting the first sample of calibration data as 
1 . For each calibration sample considered, the range of values it 
assumed is indicated by a string of asterisks, which show the limits of 
the mean sample value plus and minus the standard deviation. These 
plots show a stable but non-linear response characteristic over the 
region of the ordinate where the ordinate scale is linear. (Note that 
on both ends of the ordinate, the scale of the plot is compressed with 
respect to the central region. This was done in order to display cali- 
bration sample values over the full range of the wedge without producing 
an excessively large plot.) Over the low value range of the calibration 
data, and for all detectors, the calibration sample values are stable to 
within a standard deviation of plus or minus one count, and there is 
evident no characteristic which suggests that the calibration data is 
inadequate or inferior for any particular detector. 
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Figure 4.1-1. Plot of Decompressed Calibration Wedge Data Values and their 
Standard Deviations for Detector 1, Band 1 
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Figure 4.1-2. Plot of Decompressed Calibration Wedge Data Values and their 
Standard Deviations for Detector 2, Band 1 
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Figure 4.1-3. Plot of Decompressed Calibration Wedge Data Values and their 
Standard Deviations for Detector 3, Band 1 
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Figure 4.1-4. Plot of Decompressed Calibration Wedge Data Values and their 
Standard Deviations for Detector 4, Band 1 
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Figure 4.1-5. Plot of Decompressed Calibration Wedge Data Values and their 
Standard Deviations for Detector 5, Band 1 
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Figure 4.1—6. Plot of Decompressed Calibration Wedge Data Values and their 
Standard Deviations for Detector 6, Band 1 



4.2 GAIN AND BIAS COEFFICIENT TIME SERIES 


The MSS calibration procedure was applied to the standard six points 

from the calibration wedge data, and the resulting time series of gain 

and bias coefficients was plotted. These plots are presented in Figures 4.2-1 

through 4.2-12. In these plots, the scale of the abscissa is the same 

among all the individual gain plots and among all the individual bias 

plots , so that visual comparison among corresponding plots for different 

detectors is not deceptive. 

For the gain coefficient plots, the plotted values have been normalized 
into the range 0 to 1 , where 0 corresponds to an actual gain value of 
95.52, and 1 corresponds to an actual gain value of 122.64. Thus, the 
quantization of the plotting for the abscissa is an interval of 0.27 in 
the value of the gain coefficient. 

For the bias coefficient plots, a similar normalization has been employed 
In this case, the minimum plot value of zero corresponds to a bias value 
of -47.670563, and the maximum plot value of 1 corresponds to a bias 
coefficient value of 52.329422. Thus, the quantization of the plotting 
on the abscissa is an interval of 1.00. 

To provide an intuitive appreciation for the meaning of the fluctuations 
observed in these plots, the way in which the gain and bias coefficients 
are employed in radiometric calibration of MSS data must be considered 
The equation which is used to calibrate the data is: 



where the symbols are defined in Appendix B. This can be 
rewritten in the form: 

V = GV + B 
c 
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where 


and 


B 


■ fe) {^) 

j(sH) ^ “mu. 


+ r 


For the radiometric calibration implemented in the MDP for MSS band 1 
(0.5 to 0.6 micrometers), the nominal values for the constants in these 
expressions are: 


K 

AR q 
P = 1 





'^Max 


127 


Using these values for the constants in the calibration function, the 
expressions for G and B become: 


G = 


127 

b 

s 



s 


Since we are using unsmoothed values for the gain and bias values in the 
present investigation, the subscript 's' will be dropped, yielding as 
the expression for the calibration function: 


V 

c 


127 


CV-a) 


b 


4-10 



The effect of the observed variations in the gain, and bias coefficient 
values in the plots can be demonstrated by means of a Taylor series 
expansion of the calibration function about nominal coefficient values. 
The variation in the calibrated sample value is given by: 



f (V-a) 


Lb J 

* b=b 

n 


a=a 

o 


Aa 

b=b 

o 


127 

.2 


(Ab)V + 


127 r 

t 


a Ab - 
o 


b^a 


= (A6)v + Ab 


where a and b are typical values for the MSS calibration coefficients 
o 0 

and 

AG = - ^ (Ab) 

b 

o 

and 

Ab = ^ fa Ab - b Aa 

b2 L o 

o 

Table 4. 2-1 presents a compilation of the maximum and minimum gain 
coefficient values observed in the gain coefficient plots, together with 
the corresponding bias coefficient values, typical coefficient values 
obtained by averaging the values from the plots, and deviations from the 
typical values obtained by taking one-half the difference of the values 
obtained from the ploCs. The effect of these typical values and the 
observed deviations is shown in Table 4.2-2, which presents the values 
of the coefficients of the calibration function and series expansion of 
the calibration function, taking into account the observed inverse 
relation between variations in the gain coefficient and variations in 
the bias coefficient. 
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Note that the variation in the additive constant is limited to approxi- 
mately + 1, and this for the unsmoothed data. Thus, the variations in 
calibrated sample values attributable to this term in the calibration 
function are at most of this magnitude. The amount of variation in 
calibrated sample values attributable to the multiplicative factor 
depends on the value of the uncalibrated sample, but for a typical 
uncalibrated sample value of 20, is of the same magnitude as that due to 
the additive term. 

Pairwise comparison of the gain plots (or bias plots) for the detectors 
(e.g., detector 1 against detector 5) reveals no apparent correlation 
in the variation of the coefficient values along the time series. If 
striping, which on close inspection exhibits a six-line periodicity and 
18 therefore necessarily, even if non-causally, detector related^ were 
attributable to a systematic, detector-related variation in the cali- 
bration function, one would expect to observe a consistency in these 
variations across detectors (e.g., the bias value for detector 1 varying 
in synchronization with that for detector 5) . 

Thus, while there exist fluctuations in the calibration coefficient 
values which are of sufficient magnitude to produce a striping effect, 
these variations do not have the interdetector correlation which would 
be necessary to identify inadequate radiometric calibration processing 
procedures as the cause of the striping of MSS images. 
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Table 4.2-1. Summary of Variations in Gain and Bias Coefficients 
Obtained with Mss Calibration Procedure 


Detector 

Minimum 

Gaxn 

Bias 

Record 

Maximum 

Gain 

Bias 

Record 

Gain 

Bias 

H 

Deviation 

from 

Typical 

Ab 

Typical 

a 

o 

Deviation 

fron 

Typical 

Aa 

1 

95.51874 

1.72894 

123 

102.64764 

0.76947 

115 

99.08 

+3 56 

1.25 

+0.48 

2 

104.12521 

3.38972 

8 

114.47604 

2.86726 

71 

109.30 

+5.18 

3.13 

^,26 

3 

104.42088 

2.57248 



0.98093 

47 

108.50 

+4.08 

' 1.78 

+0.80 

4 

112.91699 

3.32911 

127 

122.64281 

1.55055 

49 

117.78 

+4.86 

2.44 

+0.89 

5 

97.20721 

4.04221 

146 

105.33539 

2.26438 

205 

101.27 

+4 06 

3.15 

+0.89 

6 

106.34644 

2.89016 

40 

114.60645 

1.52419 

197 

110.48 

+4.13 

2.21 

+0.68 
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Table 4.2-2. Variations of Calibration Function Coefficients 



MSS Calibration Coefficient Values 

Calibration Function Coefficients 


Gain 

Bias 

Gain 

Bias 

Detector 

Typical 

Value 

b 

0 

Deviation 

from 

Typical 

Ab 

Typical 

Value 

a 

0 

Deviation 

from 

Typical 

Aa 

Typical 

Value 

G 

Deviation 

from 

Typical 

aG 

Typical 

Value 

B 

Deviation 

from 

Typical 

AB 

i 

99.08 

+3.56 

1.25 

+0.48 

1.282 

+4.60x10“^ 

-1.602 

+0.673 

2 

109.30 

±5.18 

3.13 

+0.26 

1.162 

+5.51x10"^ 

-3,637 

^.474 

3 

108.50 

+4.08 

1.78 

«.80 

1.170 

+4.40x10“^ 

-2.084 

+1.014 

4 

117-78 

+4.86 

2.44 

+0.89 

1.078 

+4.50x10“^ 

-2.631 

+1.068 

5 

101.27 

+4.06 

3.15 

+0.89 

1.254 

— -2 
+5.03x10 ^ 

-3.950 

+1.274 

6 

110.48 

+4.13 

2.21 

^.68 

1.150 

+4.30x10"^ 

-2.540 

4G.877 


G => + 


127 


aG =» -Ab 


o 
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b 0 
o 


AB = 


127 

.2 


[ a Ab - b Aa I 
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4.3 INTERDETECTOR CORRELATION PLOTS 


Interdetector correlation plots of the means and root-mean-squared 
variation for 50 calibration wedge values were generated in order to 
display the relative radiometric characteristics of the six detectors 
and to evaluate the adequacy of the present linear MSS calibration 
procedure, as well as that of the various linear striping removal 
procedures. These plots use sample position within the calibration wedge 
to identify corresponding detector response samples. To the extent that 
the calibration radiances presented to each detector at corresponding 
sample positions are the same, the plots present the relationship which 
exists between the responses of the detector pairs. These plots are 
presented as Figures 4.3-1 through 4.3-5. A list of the sample positions 
in the calibration wedge data used to provide data for the plots is 
given in Table 4.3-1. Position 1 is the first sample of the calibration 
wedge data. 

These plots are essentially two-dimensional histograms where the abscissa 
scale consists of radiometric values for the reference detector (detector 1), 
and the ordinate scale consists of radiometric values for detector n (n 
= 2, 3, 4, 5, 6). In constructing this plot, the mean value of calibration 
wedge sample k of detector 1, mean value of calibration 

sample k of detector n, i^j^, plus and minus the corresponding standard 
deviations about these means, and are used to specify a region 

of the plot array, and the counts in the array elements within this 
region are incremented by 1. The numbers in the interdetector correlation 
plots are counts thus accumulated. This sort of plot provides a graphic 
representation of the functional relationship between the radiometric 
response characteristics of the separate MSS detectors (e.g., for two 
detectors viewing a set of equal radiances, if the relationship between 
the radiometric responses of the two detectors is linear, the plot will 
be a straight line; if the responses of the two detectors are the same, 
the line will have unity slope and zero intercept) • 
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Table 4.3-1. Positions in Calibration Wedge Data of Samples 
Represented In Interdetector Correlation Plots 
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All five of these plots show a relationship between detector responses 
which is essentially linear over the range from 0 to 80 along the reference 
detector axis, with a minor departure from this linearity from 80 to 110 
for all detectors. Detectors 4 and 5 exhibit the most severe departures 
from a linear relationship, this amounting to 3 counts displacement from 
a linear trend at a reference detector value of 110. 


In order to provide a quantitative comparison of the functional relation- 
ship among the six detectors the slopes, p^ = Ai^/Ai, and intercepts, 
q^, for each of the correlation plots over the linear region, from 0 to 
80 with respect to detector 1, were calculated for straight lines manually 
laid on the plots. These slopes are presented in Table 4.3-2. Also 
presented in this table are corresponding slope and intercept values 
calculated from the typical gain and bias values (B and G) given in 
Table 4.2-2, which were obtained from the gain plots of the standard MSS 
calibration coefficients. To calculate these latter slope and intercept 
values, the linear calibration functions for detector 1 and detector n 


"i = *1 b ‘i 


and 


cr = g i + b 
n “n n n 


where 


cr IS a calibrated sample value 

i^ is an uncalibrated sample value 

gj is the gain coefficient G for detector i 

b. is the bias coefficient B for detector i 
1 


are used. Note that the calibration processing for a detector is performed 
so that, for a given incident detector excitation, r, the calibrated 
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Table 4.3-2. Functional Interdetector Relationship 
for MSS Calibration Wedge Data 
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i = p i + q 
n n n n 
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sample value is independent of the detector. Therefore, it should be 
the case that: 


= §1 if (r) + bj = cr 

and cr = g i (r) + b = O" 

n ®n n n 

Thus independent of the incident detector excitation and calibrated 
sample value 


g i + b = g.i- + b 
“n n n 11 1 


or 


. H - '’n 

L * — i + 

n g n g 

n 


= p: i + q: 


n n 


n 


Using the typical values established for the calibration function gain 
and bias for as and b^, the slope and intercept of the interdetector 
relationship can be calculated. 


Figures 4.3-6 through 4,3-10 are plots of the lines whose slopes and 
intercepts are given in Table 4.2-2, comparing the interdetector relation- 
ship obtained from the calibration data with that obtained from the MSS 
calibration procedure. In these plots, the solid line presents the 
interdetector relationship established on the basis of the correlation 
plots, of calibration data, while the dashed line presents the inter- 
detector relationship which MSS calibration processing anticipates. As 
can be seen from these plots , there is substantial agreement between the 
two different interdetector relationships for detectors 3, 4, and 6, 
while for detectors 2 and 5 the MSS calibration relationship anticipates 
a different relationship from that expected on the basis of the calibra- 
tion data. 


It is tempting to identify the discrepancy observed for detectors 2 and 
5 as the cause of the observed striping. However, visual inspection of 
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the data in the smaller, low frequency region employed in the power 
spectrum analysis reveals that the striping evident in the data set 
produced by using the MSS calibration procedure on the calibration data 
for a full scene of uncalibrated data (type b2 radiometric processing) 
consists of an apparent discrepancy in the radiometric data values of 
detectors 1 and 3 (approximately level 19) with respect to detectors 2, 

4, 5, and 6 (approximately level 21). This observation is inconsistent 
with the hypothesis that a faulty calibration procedure is the cause of 
striping in data from detectors 2 and 5. 

The foregoing analysis indicates that while the MSS calibration proce- 
dure does not produce a completely accurate characterization of the 
relationships among the MSS detectors, as represented by the calibration 
data reported for the detectors, the phenomenon of striping cannot be 
attributed to such inaccuracies. It must also be noted that the validity 
of this analysis depends critically on the assumption that the calibra- 
tion wedge data used to establish the interdetector relationships repre- 
sents the response of the detector pairs to the same set of radiances. 

If the calibration lamps are not uniform, or the neutral density filter 
used to generate the wedge is not uniform across the detectors , then 
the interdetector plots would need to be compensated for any such non- 
uniformity before valid conclusions could be drawn. 
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Figure 4.3~1. Interdetector Correlation Plot for Corresponding 
Decompressed Calibration Wedge Samples and their 
Standard Deviations, Detector 2 versus Detector 1, 
Band 1 
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Figure 4.3-2. 


Interdetector Correlation Plot for Corresponding 
Decompressed Calibration Wedge Samples and their 
Standard Deviations, Detector 3 versus Detector 1, 
Band 1 
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Figure 4.3-3. Interdetector Correlation Plot for Corresponding 
Decompressed Calibration Wedge Samples and their 
Standard Deviations, Detector 4 versus Detector 1, 
Band 1 
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Figure 4.3-4. 


Inter detector Correlation Plot for Corresponding 
Decompressed Calibration Wedge Samples and their 
Standard Deviations, Detector 5 versus Detector 1, 
Band 1 
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Figure 4.3-5. Interdetector Correlation Plot for Corresponding 
Decompressed Calibration Wedge Samples and their 
Standard Deviations, Detector 6 versus Detector 1, 
Band I 
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Detector 2 Radiometric Response 
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Figure 4.3-6. Interdetector Functional Relationships for Detector 1 
and Detector 2, Band I 
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Detector 3 Radiometric Response 
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Figure 4.3-7. Interdetector Functional Relationships for Detector 1 
and Detector 3, Band 1 
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Detector 4 Radiometric Response 
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Detector 5 Radiometric Response 



Detector 1 Radiometric Response 


Figure 4.3-9. Interdetector Functional Relationships for Detector 1 
and Detector 5 , Band 1 
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Detector 6 Radiometric Response 
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A. 4 SUMMARY OF CALIBRATION DATA EVALUATION 


The analyses of the MSS calibration wedge data support the following 
4.4 SUMMARY OF CALIBRATION DATA EVALUATION 

The analyses of the MSS calibration wedge data support the following 
conclusions : 

a. Calibration wedge sample values are stable to within + 1 
count (nns) for a given sample position in the reported data. 

b. Variation in the bias adjustment in radiometric calibration of 
MSS detectors is at most + 1 count for unsmoothed calibration 
constants . 

c. Variation in the gain adjustment in radiometric calibration of 
MSS data is approximately + 1 part in 20 for unsmoothed calibra- 
tion constants. 

d. There is no evidence of any correlation in the fluctuation of 
the calibration constants which would indicate this as the 
cause of striping. 

e. The relationship between the response characteristics of any 
pair of MSS detectors, as a function of calibration wedge 
sample position, is linear over a wide dynamic range. This 
implies that linear radiometric transformations should be 
adequate to remove striping by equalizing the effective 
radiometric responses of the individual detectors in a given 
spectral band, unless, improbably, this linearity results 
from the cancellation of detector non-linearities by cali- 
bration wedge non-linearities. 
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f. Correlation plots of MSS calibration data reveal some inconsis- 
tencies between the linear interdetector relationship derived 
from the calibration data and that derived from the MSS cali- 
bration procedure. However, these inconsistencies cannot be 
identified as the cause of striping. 

The analyses of calibration wedge data have failed to establish any 
causal link between MSS calibration procedures and the phenomenon of 
striping. One additional item worthy of investigation which could not 
be pursued in this study is the effect of applying the radiometric cali- 
bration developed from the MSS calibration procedure to the calibration 
wedge data itself. This study has already established that linear 
relationships exist among the radiometric responses of the MSS detectors 
as represented by the calibration data. If the MSS calibration procedure 
is adequate within the limits of the calibration data, and if the assump- 
tion is valid that, for all detectors, the same correspondence exists 
between calibration wedge sample number and radiance incident on the 
detector, then striping should be absent from the calibrated calibration 
data, and it would thus be firmly established that the cause of striping 
is not faulty calibration processing. 

By concentrating on use of the calibration wedge data, one has a set of 
detector excitations whose properties uniquely suit them to analysis of 
interdector variations. Histogram shapes can reasonably be expected to 
be consistent without making assumptions regarding the statistical 
nature of image data. Along-track power spectra for different sample 
positions in the calibration wedge would differ primarily in the zero- 
frequency term, thus facilitating the identification of detector related, 
periodic, variations in response. With calibration wedge data, the 
along-track extent of the data sample used to produce a power spectrum 
is not limited by the size of imaged features, so that a finer frequency 
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resolution can be obtained. Finally, visual evaluation of the effective- 
ness of a radiometric processing techniques xs performed for a standard 
pattern, so that this form of evaluation is not complicated by evaluator 
biases in selecting a particular section of the image on which to base 
a decision. 
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Section 5 

INTERDETECTOR RESPONSE EQUALIZATION IN THE MDP 


Analysis of the throughput reduction and microcode changes required to 
implement the selected destriping procedure in the MDP was performed. 

The results of this analysis are presented here. 

The results presented in Section 3 of this report identified the t 3 q>e a2 
radiometric processing of uncalibrated scene data as the best procedure 
for providing interdetector response equalization for MSS scene data. 

In this procedure, detector-specific histograms of the image data from 
the first 60 mirror sweeps of the input image data for a scene are 
generated and then characterized by their means and standard deviations. 
These 6 means and 6 standard deviations are then averaged to produce a 
target mean and standard deviation. For each detector, a linear gain 
and bias radiometric transformation is then determined which, when 
applied to the scene data from that detector, will transform it so the 
mean and standard deviation of the histogram of the transformed data are 
equal to the target mean and standard deviation. This section discusses 
the effect ot incorporating this radiometric processing procedure in the 
Master Data Processor (MDP) 

The radiometric calibration of MSS data which is now performed by the 
MDP IS conceptually a two-step process: 

1. Generation of the calibration function parameters, in which 
for each detector, six samples of the calibration wedge data, 
which is supplied every second MSS mirror sweep, are used to 
calculate the gain and bias parameters of a linear transforma- 
tion of radiometric values for use on the scene data for that 
detector in each of the two mirror sweeps. 

2. Application of this calibration function to the scene data. 
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Extraction of the calibration wedge samples is performed during Input 
Processing of the data stream from the HDT, and the samples are stored 
in a disk file. Generation of the calibration function parameters from 
the samples in this file is performed on blocks of these sample cor- 
responding to an input data segment (i.e., the calibration samples 
corresponding to the data acquired between two successive nadir points 
in the World Reference System). These parameters are stored in a gain 
and bias file for that segment. Calibration of scene data is subsequently 
performed as an integral part of the resampling process, using parameters 
from these segment files which are identified by the Resampling Applica- 
tion Controller. 

The type a2 radiometric processing does not employ any of the calibration 
wedge data which is available from the MSS. It, instead requires the 
computation of the mean and standard deviation of detector-specific 
histograms, as outlined above and detailed in Appendix A. Using a 
modification of the notation of Appendix A, the type a2 processing 
requires computation of the values: 


K J 



where *^ijk scene data value of sample j from detector i for 

MSS mirror sweep k 
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J is the number of samples in the line 

I is the number of detectors used to collect scene data in 

one MSS band (6, if the thermal band is excluded) 

K is the number of MSS mirror sweeps contained in the 

subimage employed for equalization (60 in the case of 
type a2 processing) . 


The quantities 


J 


E 


u 


ijk 


j=l ^ 


are already available to MDP, since they are contained in“the 23 quality 
data pixels at the beginning of each major frame on the HDT (HDT-FM) , 
and thus are available for each detector on a line basis. The MDP-MSS 
system could be modified to strip off these sums instead of the six 
calibration sample pixel values. 

In the present MDP-MSS system, generation of the calibration function 
parameters is an autonomous task which is performed on an input segment 
basis, and which has no information available relating the specific 
input segment to an output frame. Because the type a2 radiometric 
processing requires interdetector equalization over the first 60 mirror 
sweeps of a frame, replacement of the present calibration parameter 
generation code with code to provide this radiometric equalization is 
not possible without considerable redesign of the MDP-MSS system. 
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However, if the equalization processing is performed on the scene data 
from the last 60 mirror sweeps in an input segment, rather than on the 
first 60 sweeps of the input data for an output frame, such code replace- 
ment is possible. This modification of the type a2 radiometric processing 
implies the reasonable assumption that any 60 sweep sample of the scene 
data for a frame is statistically equivalent to the first 60 sweeps. 

The radiometric equalization transformations obtained from this processing 
(one gain and bias set for each detector for each input data segment) 
would then be used for the output data frame containing the nadir reached 
at the end of the input data segment. It should be noted that this pro- 
cedure, which requires only one gain and bias set per detector per 
output frame, may decrease the size of the gain and bias file, which at 
present is required to contain a gain and bias set for each detector for 
each sweep in a swath of MSS data. 

In addition to the modifications discussed above, the present MDP-MSS 
system would have to be modified so that the Resampling Application 
Controller would cause the same gain and bias set per detector to be 
applied over an entire output data frame. 

Since the changes required to perform this form of radiometric processing 
for destriping on the MDP-MSS system can be accomplished as described 
above, a procedure which does not require any major revisions to the 
system data flow, the system throughput can be expected to be unchanged. 



Appendix A 

IBM DESTRIPING PROCEDURES 


A.l SWEEP HISTOGRAM EQUALIZATION 


An algorithm of the following type satisfies the requirement for adjusting 
detector gains and offsets between detectors to minimize striping. 


a. Compute a histogram for each detector. 

b. Adjust the gain and offset of each detector so that the mean 
and standard deviation of the histogram of its corrected 
output matches the average of the means and standard deviations 
of histograms produced for each detector. 


Let the average mean and standard deviation for the detector histograms 
M, S: let the mean and standard deviation for the ith detector be ra., 

, and let the output of the ith detector for sample j be u_ . The 
following equation then transforms u^^. into a variable w_ having the 
mean and standard deviation of the ensemble: 


"ij • “ '“ij - “i> 


The corrected output for the ith detector is w^^. . The gain and offset 
for the ith detector then are: 



The quantities m^, s^, M, and S should be evaluated over a span of data 
such that the ground area scanned by each detector is statistically 
equivalent. During a single sweep, the set of inputs to each detector 
is assumed to meet this criterion. 


Let u . be the jth sample from the ith detector. 


Compute for i=l,2,3,4,5, 


m 


i 


1 ^ 

n 




(n = number of samples in the line) 



These computations can be implemented as follows. While a scan line is 
being processed, the sum and the sum of squares are accumulated and used 
to calculate and values to be used on a following sweep. In the 
same pass the corrected values are calculated using the following 
equation: 


w. . = B. + G. u . 

ij 1 1 ij 

The values of and used are from the data of the sweep previous to 
the sweep being corrected. This technique avoids having to make two 
passes of each scan line: one to compute the and values and 

another to use those and G^ values to modify the line. 



A. 2 SUBIMAGE HISTOGRAM EQUALIZATION 


In the procedure for sweep histogram equalization described above, the 
assumption is made that the data provided by any detector for one mirror 
sweep IS statistically equivalent to that provided by the other detectors 
in the same sweep. An alternative procedure, which has the advantage of 
increasing the data on which this statistical equivalence is imposed, is 
to base the interdetector equalization on the image data from some 
significant fraction of the entire image* This image fraction can be 
anywhere from one sweep to the entire image. 

In this alternative equalization procedure, the gain and bias radiometric 
constants for each detector are calculated and outlined above, but based 
on histograms derived from all the image data contained in a specified 
subimage. The radiometric constants are then applied to the entire image 
data set. 
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Appendix B 

MSS CALIBRATION PROCEDURE 


The radiometric calibration formula employed for the four MSS bands 
available on Landsat 1 and 2 is: 



where = calibrated pixel value; 

V = input (decompressed) pixel value; 

V^ax ” maximum pixel value (127 for three bands and 63 for the fourth 
band) ; 

Afi = R - R . , depending on each band; 
max min’ f ^ 

K, p, q, and r are constants (expected to vary infrequently); 

= smoothed offset (computed per scan line); 
bg = smoothed gain or slope (computed per scan line). 

This calibration formula assumes that the data either was taken in the 
linear mode or has been decompressed. 

The offset and gain coefficients are determined once per scan line as 
follows : 

a = £ V b = XI D^ 

i i 

where i runs from 1 to 6. 
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The are the linear or linerized detector calibration samples associated 
with the scan line; C^, are regression coefficients, which may differ 
for each detector and band, but remain constant during normal sensor 
performance. 

Smoothed offset and gain coefficients, a and b , for each detector are 

s s 

calculated for every mirror sweep, n, as follows: 

b^(n) ® bg(n-l) + W(n) 

a^(n) “ a^(n-l) + W(n) ja(n) - a^(n-l)j 


where 

n = sequential number of mirror sweep; 

b(n), a(n) = latest value of b and a determined; 

b (n) , a (n) = new smoothed value of b and a; 
s s 

W(n) = 1/n+l for n=l to 15, 

W(n) = 1/16 for n ^ 16. 

The and are predetermined on the basis of prelaunch radiance test 
data. 
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Appendix C 
PLATES 


Plate 1 
Plate 2 

Plate 3 

Plate 4 

Plate 5 

Plate 6 

Plate 7 


- MSS Scence 2183-16433, MSS Band 1, Uncalibrated Data 

- MSS Scence 2183-16433, MSS Band 1, Uncalibrated Data after 
Type al Radiometric Processing 

- MSS Scene 2183-16433, MSS Band 1, Uncalibrated Data after 
Type a2 Radiometric Processing 

- MSS Scene 2183-16433, MSS Band 1, Uncalibrated Data after 
Type a3 Radiometer Processing 

- MSS Scene 2183-16433, MSS Band 1, Uncalibrated Data after 
Type bl Radiometric Processing 

- MSS Scene 2183-16433, MSS Band 1, Uncalibrated Data after 
Type b2 Radiometric Processing 

- Scene Overlay Showing Specific Regions Employed in Evaluations 


(These plates are contained in the envelope at the end of this report) 




Plate 1 - MSS Scene 2183-16433, MSS Band 1, 
Uncalibrated Data 





Plate 2 - MSS Scene 2183-16433, MSS Band 1, 
Uncalibrated Data after Type al 
Radiometric Processing 








Plate 3 - MSS Scene 2183-16433, MSS Band 1, Uncalibrated 
Data after Type a2 Radiometric Processing 

















Plate 4 - MSS Scene 2183-16433, MSS Band 1, Uncalibrated 
Data after Type a3 Radiometric Processing 






Plate 5 - MSS Scene 2183-16433, MSS Band 1, Uncalibrated 
Data after Type bl Radiometric Processing 

















Plate 6 - MSS Scene 2183-16433, MSS Band 1, Uncalibrated 
Data after Type b2 Radiometric Processing 
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Plate 7 - Scene Overlay Showing Specific Regions Employed in Evaluations 



